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Abstract 

We study the general-setting quantum geometric phase acquired by a particle 
in a vibrating cavity. Solving the two-level theory with the rotating-wave 
approximation and the SU(2) method, we obtain analytic formulae that give 
excellent descriptions of the geometric phase, energy, and wavefunction of 
the resonating system. In particular, we observe a sudden 7r-jump in the 
geometric phase when the system is in resonance. We found similar behaviors 
in the geometric phase of a spin-1/2 particle in a rotating magnetic field, for 
which we developed a geometrical model to help visualize its evolution. 
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I. INTRODUCTION 



The dynamics of confined cavity fields interacting with the cavity wall is of great interest 
for the understanding of a variety of problems such as hadron bag models [1] , sonolumines- 
cence [2], cavity QED [3] and black hole radiations [4]. Most previous studies on dynamical 
cavity concentrated on scalar or photon fields [5], which despite the apparent simplicity, 
exhibit rich and complex dynamics. In contrast, the system of a Schrodinger particle in an 
oscillating spherical cavity has not gained as much attention. In this article, we study the 
generalized quantum geometric phase of a particle in resonance with the vibrating cavity. 
We will show that the geometric phase acquires sudden 7r-jumps when the particle makes 
transitions. 

The geometric phase of a quantum system has drawn much attention since it was first 
proposed by Berry in 1984 [6,7]. It reveals the geometric structure of the Hilbert space of 
the system [7,8], and its existence in physical systems has been verified in many experiments 
[9], including electron diffraction from a screw dislocation [10] and neutron interferometry 
[11]. The phase effects in molecular physics can lead to energy splittings and shift quantum 
numbers. The geometric phase has also been shown to be intimately connected to the physics 
of fractional statistics, the quantized Hall effect, and anomalies in gauge theory [7]. As far 
as we know, our study represents the first calculation of the geometric phase of a resonating 
system, which evolves non-adiabatically and non-cyclically. 



II. FORMALISM 

We consider an infinite cylindrical or spherical potential well with oscillating boundary 
[12]: 

where R(t) = Ro/a(t), with a(t) = [1 + esin(u;£)] _1 . The coordinates can be transformed to 
a fixed domain via y = a(t)r, and, to preserve unitarity, the wavefunction is renormalized 
through <f)(y,t) = a~^ip(r,t), where £ = 1 (3/2), for a cylindrical (spherical) well. 

Since the full Hamiltonian H(t) commutes with L 2 and L, the wavefunction can be 
factorized: (fi(y,t) = Y(y,t)Q, where 9 depends only on the angular variables. Inside the 
cavity, the radial wavefunction satisfies 



dY 
~dt 



iho? ( d 2 | ndd_ _ m d \ + R ( d_ + ^ N 
2/i \dy 2 y dy y 2 I R \ dy 



Y, (2) 



where /i is the particle mass, and n d — 1, 2 and = m 2 , 1(1 + 1) for cylindrical and 
spherical wells respectively. In this paper, we only consider the rrid = sector. 

The wavefunction described by Eq. 2 evolves in time and acquires a time-dependent 
phase, which in general consists of a dynamical phase and a geometric one [6]. When the 
dynamics is adiabatic and/or the evolution is cyclic, the geometric phase, or Berry's phase, 
has been studied for many systems. Since we are interested in the geometric phase of a non- 
adiabatic, non-cyclic system, we have to resort to a generalized method. Following Ref. [13], 
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we first remove the dynamical phase from the wavefunction of the system: \(f>(t)) — > \(f*(t)) = 



-iO{t) 



\4>{t)) , the dynamical phase is 



0(t) = -\f Q E(t')dt' , 



(3) 



where E(t') = Re(<P{t')\H(t')\<p(t')) . 

The general setting geometric phase, or the Pancharatnam phase f3 is defined as the 
relative phase between the state at time t\ and that at time t, both with the dynamical 
phase removed. It can be obtained from the inner product of these two states 



P{h,t) = -Hog 



(4) 



Our main goal in this paper is to study this geometric phase for a simple but nontrivial 
dynamical system. 



III. RESULTS 

We have solved Eq. 2 numerically and checked that the solution is stable and converges 
very well. In Fig. 1, we show the maximum energy of the particle as a function of the 
driving frequency u, both having been scaled to be dimensionless: E m3iX = fiRlE maiX /h 2 , 
and lo = iiRqUj/H. The particle is initially in the ground state, and it is in resonance at 
specific values of Co. In each period of cavity vibration r = 2ir/Co, the geometric phase 
acquires some changes, (3\(ti) = £i +r), and it exhibits oscillations at these resonances. 
We show three examples in Fig. 2 for Co = 66.632, 17.278 and 22.21227, where the geometric 
phase acquires periodic changes of 27r, 4n, 6n respectively. All the resonances shown in 
Fig. 1 are associated with oscillations of Pi with amplitudes of 2Nir, N — 1, 2, 3. 

To understand the resonances and the associated geometric phases, we limit ourselves 
to the parameter regime where the Hamiltonian can be truncated to a two-level system. 
Specifically, if e << 1, and the initial state is the k th unperturbed eigenstate \4>k{y)), with 
eigenenergy Ek, then when the driving frequency corresponds to the energy difference be- 
tween the initial state and another unperturbed eigenstate |0„(y)), with eigenenergy E n , i.e., 
lo = uo n k = (E n — Ek)/h, the particle is expected to behave as in a two- level system. Then 
the problem simplifies considerably, and we have obtained its solution with two analytic 
approaches: the SU(2) method and the Rotating- wave Approximation (RWA). 

Following Cheng et al. [14] we first expand the time- dependent Hamiltonian H(i) of a 
two- level system in the identity operator J, the raising and lowering operators a±, and Pauli 
spin matrix <r 3 

H(t) = i [f (t)I + h(t)a + + / 2 (t)<7 3 + h(t)a.] , (5) 

where fj, j = 0,1,2,3, are in general complex functions of time. The evolution operator 
can then be written as 

U(t, 0) = b(t)eRWW ( b ~ 2{t) + g3 9 { f )93(t) ^ ) , (6) 
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where b(t) = exp(— g2(t)), and gi satisfies 



g\ = h + 2/ 2 #i - f 3 gf , 
92 = h ~ h9i , 
gs = fsb~ 2 , 



(7) 
(8) 
(9) 



with the initial conditions gj(0) = for j =1, 2, or 3. 

Suppose that the two levels in which the system oscillates are \4>k{y)) and \(f> n {y)), so 
that the wavefunction is \(f)(y,t)) = ak(t)\<fik(y)) + On(^) 1 0n (?/))• Furthermore, if the initial 
conditions are 0^(0) = 1 and a n (0) = 0, then 



afc W ^ =m JlHt')dt' ( b 2 + gig 3 \ 
a n (t) J V 9s J ' 



(10) 



The SU(2) method is exact and reduces the problem to solving the ODE's for gj. 

The spirit of the RWA is to retain only those terms in the Hamiltonian that correspond to 
the resonance frequency. We first separate out the fast phase factors from the wavefunction: 

\<P(y,t)) = c^e-^lMy)) + c B (t)e-*" (t) |&(v)>, (11) 
with piif) = Ei Jq a 2 (t')dt' jfi. Substituting Eq. 11 in the Schrodinger equation, we have 



Cfc 



- W *(t) \ ( c k 
W{t) [c n 



(12) 



where 



^ffexp 



f a 2 (t')dt' 
Jo 



and 



Vnk 



4>n 



(13) 



(14) 



is a constant that depends only on the states involved. 

As e is small, we can expand W(t) as a series of e. Retaining all terms up to third-order, 
we have: 



W(t) 



ee 



iu)„, k t 



iu n k{cos 2ut + 1) — — sin 2ut 



+ 
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uj cos cut + e 
{uj 2 — QuJni^j cos cot — (uj 2 + iJ 1 ^ cos 3uit 



— i-(j n k<jj(smut + sin3a;t) 
2 



(15) 



Note that W(t) consists of oscillatory terms with various frequencies depending on uj. In 
the spirit of RWA, we keep only terms with the lowest frequency for each uj, the rationale 
being that and c n vary slowly in time and the contributions to W from high frequency 
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terms cancel on average over such a long time scale. It is clear then from Eq. 15 that W 
is largest if 00 = tv nk /N, with N an integer, because of the emergence of zero frequency 
terms. For these driving frequencies, W is large and effective in inducing transitions, and we 
have resonances. At or close to a resonance, Eq. 15 simplifies tremendously, and closed-form 
solutions for Eq. 12 can be obtained easily. 

For example, when Sou = u — uj nk ~ 0, W(t) ~ ecuexp (—iSut/2) /2, and we have Rabi 
oscillations [15]: 

Ck (t) = e l5 ^ 2 [cos xt - iScu sin X t/(2 X )] , (16) 
c n (t) = Te~ l5 ^ 2 sm X t/ X , (17) 



where we took Cfc (0) = 1, c n (0) = 0, x = yT 2 + 5io 2 /4, and r = Ti = rj nk eu;/2. The 
maximum value of |c n | 2 , |c n |^ ax (a;) = 1/[1 + (5uj/2T) 2 ], is a Lorentzian with a FWHM of 4r. 

The solutions for higher N resonances differ only in the widths r = = e N rj nk ^M /2N , 
where 



In 



(uj nk + 5uj) N = 1, 

(3cu nk - 5uj) /4 N = 2, (i 8 ) 

7 U>l k -17LL) nk 8u)+2Sul 2 AT _ O 

24(LU nk +Sw) iV - J > 



with 5c<j = A^cj — uo nk . The FWHM of the resonances AT /N oc are narrower for larger TV, 
and the oscillation periods T oc e"^ are longer. 
The energy of the system is given in RWA by 

E(t) w a 2 {t) [E k cos 2 + (A + E n ) sin 2 xt] , (19) 

where A = [-Efe<5a; 2 + 4i? n (r 2 — x 2 )] /4x 2 • The energy of the particle therefore exhibits fast 
but small oscillations due to the a 2 factor, modulated by a large but slow oscillation of 
period ir/x- It follows that 

(1 — e) 2 E m ^ — E k 1 2 

MLxM , (20) 



E n -E k 1 + (c^/2r) 2 

which is again a universal Lorentzian for all resonances. 

The RWA results are virtually identical with those of the SU(2) method. Fig. 1 can 
be understood completely in terms of overlapping series of the iv nk /N resonances, and the 
positions of the peaks are in almost exact agreement with the RWA or SU(2) predictions. 
The resonance line shapes are well fitted by the Lorentzian Eq. 20, though the widths 
are underpredicted for the N = 2, 3 resonances, as shown in Table 1. For example, the 
(N — 2, n — 4) resonance has a much larger width than predicted by RWA and SU(2), which 
is mainly due to the involvement of other states during transition, such as a (1 — > 3 — > 2) 
process, which is second order in perturbation theory and so affects the N > 1 transitions 
more severely. Occasionally, resonances at similar frequencies may overlap and lead to 
broadened widths. 

At a resonance, the Rabi oscillations of the particle as predicted in RWA (Eq. 19) can 
be seen explicitly in Fig. 2, where the energy of the particle vs. time is shown for n = 4, 3, 4 
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and N = 1,2,3. The RWA also predicts that at exact resonances, x — F, and therefore 
the oscillation periods should be inversely proportional to the widths of the resonances, 
T = Tf/x = tt/T. As shown in the last column of Table 1, this is clearly borne out in 
the numerical data, where TY/tx is listed and is found to be closed to 1 for all resonances, 
deviating by less than 10% for even those with strong mixing. 

Using the RWA and two-level approximation, we can calculate the dynamical phase 
easily. Removing the dynamical phase with the help of Eq. 19, we have 

m))«^^ w ( u ^L,), (2D 

where oj' nk = A/fi + oj nk . For oj = uj nk /N resonances, if we choose t — t-y to be an integral 
multiple, q, of the cavity oscillation period r, we get 

(3(h, h + qT ) = ^[x- sin X cos(2 X ti + x)} , (22) 

where x = X1 r - I n particular, if t\ = 0, 

for (2m - \) < t/T < (2m + ±), 
± n, for (2m + \) < t/T < (2m + §), 



Pom = m v) = \a^:_z_ )Z :ic. !L l ):: ic (23) 



where Q(t) = u n kt — ^sin2xt , with m = 0,1,2,.... There are sudden approximate n 
jumps in (3 at t/T = (2m + l)/2, as shown for example in Fig. 3 for Co = 12.344. Since 
XT « 1, the phase change in each cycle (q = 1) is 

A(ti) =(3{t 1 M J rr) ^u' nk rsm 2 x(t 1 + r/2). (24) 

At an exact A^ th resonance, 5lu = and A = 0, and so 

/3i(f i) w 2iV7r sin 2 + Nir/Lu nk )] . (25) 

Therefore, the geometric phase oscillates with an amplitude of 2Ntt and period of T. Both 
of these RWA predictions are in excellent agreement with the numerical data, as shown 
in Fig. 2. Note that the n -jumps and the functional forms of the geometric phases are 
independent of e and h, as long as they are nonzero. 



IV. AN ELECTRON IN A MAGNETIC FIELD 

To gain more insight into the geometric phase for a two-level system, we have studied a 
simple model of a magnetic field rotating around a spin-1/2 particle. Suppose an electron 
of charge — e and mass m is placed at the origin, in the presence of a magnetic field 

B(£) = -B [sin a cos(ut)i + sin a sin(ut)j + cos ak] , (26) 

which has a constant magnitude B but its direction sweeps out a cone with an opening 
angle a,0<o;<7r, at a constant angular speed oj. The Hamiltonian of the system is given 
by 
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H(t) — — B • S = — — ^ ( ™. a e ^ W ) , (27) 
w m 2 \ e l0Jt sin a - cos a J ' v 7 

where u;i = —eB Q /m and S is the spin matrix. 

The system can be solved analytically [16]. The instantaneous eigenspinors of H(t) with 
eigenenergies E + = —E_ = —huoi/2 respectively are 

i*w> = ( e ^£/2)) < 28 > 

and 

v-w-iJttgm)- (29) 

Suppose the electron spin is initially parallel to B(0), and we consider the case when 
the particle makes a transition to spin down along the instantaneous direction of B. This 
happens with unit probability if uj — —uji/ cos a, provided that cos a ^ 0. The state vector 
at any time is then [16] 

\m) = e^ t/2 jcos (y) \Mt))+i** (y) l^-(0>} , (30) 
where A = u sin a . 

The Pancharatnam phases comparing the initial state with the state at time t for different 
values of a are shown in Fig. (4-6). The sudden 7r-jump and a two-period oscillation can be 
seen in the case when a — 1/p, p being an even integer, such as the case in Fig. (4). When 
p is an odd integer, the Pancharatnam phase performs a single-period oscillation with no 
7r-jump, as shown in Fig. (5). It has neither single-period nor two-period oscillation when 
other values of a are used, Fig. (6). If t = qr, with q an integer and r = 2tv/uj, 

f for(2m-i)<gr/T<(2m+i), 
m ) ~ \ -2<gli ± tt, for (2m + !)< qr/T < (2m + |), {6l) 

which has the same form as Eq. 23, where m is an integer, T = 2it/\, and 

Q(t) = sin At + utj . (32) 

For the case t 1 — (q — l)r and t = qr, 

a (j \ _ f ~\ Wi +r)- fi(ti)] , for < sin a < 1/2, 

PlKl) 1 -i + r) - ^(ti)] ± tt , forl/2<sina<l. 1 } 



In the limit « C 1, 



w -2tt sin 2 ( ^(*! + r/2) ] , (34) 
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which has the same form as Eq. 25. 

In the cyclic limit, i.e., t = qr = T = 2it/\, 

tt(T) =luT = ujqr = 2qix , (35) 

and thus the geometric phases are fio{T) = — (q — and /3\(T) = —tt. 

Since the geometric phases in the two models - an electron in a rotating magnetic field and 
a particle in a vibrating cavity - are remarkably similar, we conjecture that the main features 
of the generalized geometric phase we calculated, especially the 7r-jumps, are universal for a 
particle in transition from one state to another. Similar features about the geometric phase 
have also been obtained in [17]. 



V. GEOMETRICAL MODEL 

Here we present a geometrical model that helps to visualize the evolution of the geometric 
phase of a two-level system. It is clear from Eq. 30 that the state vector traces out a path 
on the unit sphere defined by the angular coordinates 9 = Xt and = uot. Note that we 
have used a time-dependent basis \ip + (t)), which mark the North and South Poles 

respectively on the unit sphere. The excitation condition of our system therefore gives us 
the trajectory of the state on this unit sphere, 

9 = sin a , (36) 

which is a spiral curve. 

The solid angle subtended at the origin by the spiral curve and the geodesic, = 0, up 
to any t = qr is 

/•2qir p(j>' sin a 

n o = / sinflWd0' (37) 
jo Jo 

sin (2qn sin a) , . 

= 2g7r y -4 >-, 38 

sin a 

which coincides with Eq. 32, in the limit a < 1. Therefore, we see from Eq. 31 and Eq. 33 
that the geometric phases (3q and f3\ are simply the solid angles subtended by the spirial 
curve, 6 = et(f>, from = and = 2(q — l)n up to = 2q7i respectively, divided by —2. 
This picture is a generalization of Berry's [6] for adiabatic and cyclic evolution corresponding 
to trajectories with constant 9 = 9 , and the solid angle is simply £l = 2n(l — cos9 ). Even 
the 7r-jump can be represented in this model by a jump of the particle from one sphere to 
another if its trajectory happens to reach the South Pole, which depends on the value of the 
angle a. The double-sphere picture reminds one of the underlying SU(2) structure. 



VI. SUMMARY 



In summary, we have reported the first calculation of the quantum geometric phase of a 
physical system in resonance - that of a particle in a vibrating cylindrical or spherical cavity, 
and we have shown that it acquires sudden 7r-jumps when the particle makes transitions from 
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one state to another. We have derived analytic expressions in the RWA and SU(2) methods, 
which give excellent description of the energy, wavefunctions, and the geometric phases at 
these resonances. We found remarkably similar properties of the geometric phases for the 
simple system of an electron in a rotating magnetic field, which led us to conjecture that 
the main features of the generalized geometric phases and especially the 7r-jump we found 
are universal. We have also developed a geometrical model to help visualize these phases. 

This work is partially supported by a Chinese University Direct Grant (Project ID: 
2060093). 
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FIG. 1. The maximum energy of the particle as a function of the driving frequency, obtained 
by solving the Schrodinger equation numerically, for a cylindrical cavity with e = 0.01 (solid line) 
and a spherical cavity with e = 0.02 (dashed line). 
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FIG. 2. The scaled energy, S = (Ea~ 2 — E^)/{E n — Ef.) where E^ and E n are the unperturbed 
eigenenergies, and geometric phase acquired at each cycle, 0i/2Nir, vs. t/T, T being the oscillation 
period of the geometric phase, for the case of cylindrical cavity at u = 66.632, 17.278, and 22.21227, 
using e = 0.01. 
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FIG. 3. (3a vs. t/T at u = 12.344 for a cylindrical cavity with e = 0.01. The region around 
t/T = 0.5 is enlarged and shown in the lower panel to show the approximate 7r-jump. 
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FIG. 6. Same as Fig. 4, but for a = 5/501. 
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TABLES 



Resonances 


r 


TT/tt 


N 


n 


Numerical 


RWA 


Numerical 


1 


2 


6.17 (7.45) 


6.17 (7.40) 


1.00 (1.01) 


1 


3 


16.6 


17.3 


0.98 


1 


4 


33.2 


33.3 


1.00 


2 


2 


4.15 (4.66) 


4.63 (5.55) 


0.99 (0.96) 


2 


3 


50.0 (76.5) 


13.0 (14.8) 


1.00 (1.02) 


2 


4 


891 (532) 


25.0 (27.8) 


1.09 (1.17) 


2 


5 


136 


40.7 


0.93 


3 


4 


205 (341) 


23.6 (26.2) 


1.02 (1.07) 


3 


5 


15700 


38.5 


1.02 


3 


6 


78300 


56.8 


1.07 



TABLE I. Scaled widths (f = [iB^T /fie N 7]^) of the lowest k = 1 resonances for cylindrical 
cavity with e = 0.01. Full numerical results are compared to the RWA values. The widths multiplied 
by the Rabi oscillation period, both extracted from the numerical data, are also shown in the last 
column. Corresponding numbers for spherical cavity with e = 0.02 are shown in parentheses. 
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